Ab initio study of BaTiOs and PbTiOs surfaces in external electric fields 
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For the ferroelectric perovskite compounds BaTiOa and PbTiOa, we have studied the effects of 
external electric fields on the structural properties of the (001) surfaces. The field-induced changes 
in the surface interlayer spacings and bucklings have been calculated using a first-principles ultrasoft- 
pseudopotential approach, and the change of the polarization and the ferroelectric distortions in the 
surface layers have been obtained. The surfaces are represented by periodically repeated slabs, and 
an external dipole layer is included in the vacuum region of the supercells to control the electric 
field normal to the surfaces. The infiuence of the electrical boundary conditions on the ferroelectric 
properties of the slabs is discussed. In the case of a vanishing internal electric field, our study 
indicates that even very thin slabs can show a ferroelectric instability. 

PACS numbers; 77.55. +f, 68.35.Bs, 77.22. Ej, 77.84.Dy 
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I. INTRODUCTION 

In recent years, ferroelectric (FE) compounds based 
on the cubic perovskite structure ABO3 have attracted 
much interest because of their promising potential for a 
series of technological applications. Just to mention two 
examples from semiconductor industry, their ferroelectric 
properties can be used to build non-volatile random ac- 
cess memories, and their high permittivity makes them 
good candidates to replace Si02 in large-scale integrated 
circuits. Because of the ongoing miniaturization in semi- 
conductor devices, this may be soon required in order to 
be able to construct capacitors with sufficient capacity 
(for example in DRAM cells) and to maintain gate oxide 
layers in MOSFETs which are thick enough to prevent 
tunneling of electrons between the gate electrode and the 
channel. 

In these applications it is the aim to apply the FE ma- 
terials in very thin film geometries where the ferroelec- 
tric and dielectric properties will be strongly influenced 
by surface effects. For the understanding of the behav- 
ior of FE thin films it is therefore important to explore 
how the FE order parameter couples to the surface. In 
many experimental investigations of thin films, a degra- 
dation of the FE properties has been observed. It has 
been believed that decreasing the thickness of thin films 
suppresses ferroelectricity and eliminates it altogether at 
a nonzero critical thickness which has been estimated to 
be approximately 10 nm for PbTiOa.El However, in recent 
experiments, Tybell and coworkers have shown that for 
high-quality films of Pb(Zr2,Tii_2,)03, a stable polariza- 
tion perpendicular to thcj-surface persists down to film 
thicknesses of 40 A or less.0 This suggests that the previ- 
ously observed suppression of ferroelectricity in thin films 
is not a purely intrinsic effect caused solely by the pres- 
ence of the surface, but is related to extrinsic factors like 
perturbations of the chemical composition of the surface 



(impurities, oxygen vacancies or other defects), surface 
induced strains, or variations in electrical and mechani- 
cal boundary conditions. 

To investigate this further, we have studied slabs of 
isolated, ideal (i.e., clean and defect-free) thin films of 
BaTiOa and PbTiOa in (001) orientation. Both possible 
surface terminations (on BaO/PbO or Ti02 layers) have 
been considered. It is the aim of this study to determine 
whether thin films of this kind still show a ferroelectric 
instability, and how the ferroelectric distortions and the 
spontaneous polarization change at the surface. 

A few theoretical investigations in this direction have 
already appeared. In Refs. ra and H we have studied 
the (001) surfaces of BaTiOa and PbTiOa for the case 
of the tetragonal axis (i.e., the direction of the spon- 
taneous polarization) lying parallel to the surface. Al- 
though the surface relaxation energies are found to be 
substantial (i.e., many times larger than the bulk ferro- 
electric well depth) , it turned out that the influence of the 
surface upon the FE order parameter is only modest. For 
the Ti02-terminated surface of BaTiOa and the PbO- 
terminated surface of PbTiOa, a modest enhancement 
of the FE order was found at the surface; for the other 
two surfaces, a small reduction in the FE distortions was 
observed. However, in all cases the ground state was fer- 
roelectric, and deviations of the FE distortions from the 
bulk value were confined to the first few atomic surface 
layers. 

On the other hand, for thin PZT films the experimen- 
tal investigations of Tybell have shown a polarization 
perpendicular to the surface. In this case, theoretical 
calculations have to deal with the additional problem of 
the correct electrical boundary conditions far frppi the 
surfaces. In two previous studies on BaTiOaJj'El slabs 
of truncated bulk material with a net polarization per- 
pendicular to the surface have been used. The slabs were 
either repeated periodically, assuming periodic boundary 



conditions for the electrostatic potential, or were treated 
as isolated slabs with a vanishing external electric field, 
and in neither case were the atoms allowed to relax fully 
from their ideal bulk positions. 

We point out in this paper that these situations are 
rather artificial and do not provide useful information 
about the FE properties of the slabs. Instead, we will 
show that the appropriate electrical boundary condition 
is a vanishing electric field inside the slab. This implies 
that if the slab shows a spontaneous polarization, an ex- 
ternal electric field has to be applied to compensate the 
depolarization field caused by surface charges. A van- 
ishing internal electric field is equivalent to the so-called 
short-circuit boundary conditions, where the thin film is 
sandwiched between grounded plates of a capacitor. 

Short-circuit boundary conditions have been used by 
Ghosez and RabeEl in a microscopic effective-Hamiltonian 
study of PbTiOa thin films. They found that (001)- 
oriented films with thicknesses as low as three lattice 
constants have a perpendicularly polarized ferroelectric 
ground state with a significant enhancement of the polar- 
ization at the surface. However, in the effective Hamil- 
tonian, no information about the structural relaxation 
of the surface layers was included. With the present ab 
initio calculations, we confirm the finding of Ghosez and 
Rabe about the ferroelectric ground state, but we also 
show that the atomic relaxations at the surface signif- 
icantly modify the effective-Hamiltonian picture of the 
polarization in the surface layers. 

The paper is organized as follows. In Section II we de- 
scribe the technical details of our computational method 
and the geometry of the slabs. We also discuss how to de- 
duce the correct electrical boundary conditions and how 
to apply these in the supercell calculations. In Section 
HI we present the results of our calculations on various 
slabs exposed to external electric fields. Finally, the pa- 
per ends with a summary in Section IV. 



II. THEORETICAL DETAILS 

A. Method of calculation 

As in our previous studies on perovskite surfaces 
we have carried out self-consistent total-energy calcul 
tions within the framework of density-functional theorj 
using the Vanderbilt ultrasoft-pseudopotential scheme.l 
In the pseudopotential generation the semicore Ba 5s and 
5p, Pb 5d, and Ti 3s and 3p orbitals have been included 
as valence states (for more details on the pseudopoten- 
tials see Ref. ^l]). The electron wave functions were ex- 
panded in a plane-wave basis set including plane waves 
up to a cutoff energy of 25 Ry. As has been shown in 
previous studies, this cutoff energy is sufficient to obtain 
well converged results for these materials. A conjugate 
gradient technique was used tOjOiiinimize the Hohenberg- 
Kohn total energy functional,!!^ and the exchange and 



correlation effects were treated within the local density 
approximation (LDA) in the Perdpv-Zunger parameteri- 
zation of the Ceperley- Alder data.tHi The positions of the 
ions were determined by minipiizing the atomic forces 
using a variable-metric scheme J13 The surfaces were con- 
sidered to be fully relaxed when the forces on the ions 
were less than 0.01 eV/A. 



B. Surface and slab geometries 

BaTiOa and PbTiOg both belong to the group of II-IV 
perovskite compounds, i.e., ABO3 perovskites in which 
atoms A and B are divalent and tetravalent, respectively. 
In this case, the AO and BO2 layers are charge-neutral, 
so that both AO- and BO2 -terminated surfaces are non- 
polar. 

The surfaces have been represented by periodically re- 
peated slabs consisting of 7-9 alternately stacked lay- 
ers of AO and Ti02. Both types of surface termina- 
tions have been considered, and the slabs were sepa- 
rated by a vacuum region of two lattice constants. A 
(4,4,2) Monkhorst/Pack k-point meshlia was used for all 
Brillouin-zone integrations. As has been stated in Ref. |j, 
the results for structural properties of perovskite surfaces 
are very well converged for this choice of the supercell ge- 
ometry and k-point mesh. 

BaTiOs and PbTiOs display different sequences of 
structural phase transitions as the temperature is low- 
ered. PbTiOa undergoes a single transition from a cubic 
paraelectric to a tetragonal ferroelectric phase at 763 K, 
which is the ground state structure at T=0. BaTiOa dis- 
plays a series of three transitions from cubic paraelectric 
to tetragonal, orthorhombic, and rhombohedral ferroelec- 
tric phases at 403 K, 278 K, and 183 K, respectively. 

Because we are primarily interested in the situation 
where a material shows a polarization perpendicular to 
the surface, we focus here on the tetragonal FE phases of 
BaTiOa and PbTiOs, and we consider only the case of the 
tetragonal c axis pointing perpendicular to the surface 
(in the following referred to as the z direction). To pre- 
vent BaTiOa from adopting the true rhombohedral T—0 
structure, mirror symmetries M^, and Mj, have been im- 
posed during the relaxation of the atoms, thus mimicking 
the experimental room-temperature structure. In other 
words, only displacements of the atoms in the z direction 
(perpendicular to the surface) were allowed. This also 
prevents the polarization from rotating to become paral- 
lel to the surface. Additionally, the slab lattice constant 
in the x and y directions was set equal to the theoreti- 
cal equilibrium lattice constant a computed for the bulk 
tetragonal phase (0=3.94 A and 0=3.99 A for BaTiOa, 
and a=3.86 A and c=4.04 A for PbTiOs). 



C. Electrical boundary conditions 

For slabs with a net polarization perpendicular to the 
surface, two questions have to be addressed: (i) What are 
the appropriate boundary conditions for the electrostatic 
potential in the FE state? (ii) How will the electrostatic 
potential be modified by the periodic repetition of the 
slabs? 



1. Boundary conditions in the FE state 

Let US first consider an isolated slab with a polariza- 
tion perpendicular to the surfaces. We choose the surface 
normal n to be parallel to the z-axis, and we assume the 
charge density p{r) of the slab to be periodic in the x- 
and y-direction. The polarized slab exhibits an electric 
dipole moment 



p{z) z dz 



(1) 



parallel to the surface normal n, where p(z) is the planar 
averaged charge density 



AA 



— p(r) dxdy 



(2) 



and A is the area of the surface unit cell. The electro- 
static potential v{r) seen by the electrons can be calcu- 
lated by solving the Poisson equation 



V^w(r) = 47rep(r) 



(3) 



In addition to the microscopic quantities p{r) and v{r), 
we assume the slabs to be thick enough that also macro- 
scopic quantities like the macroscopic electric field E, the 
dielectric displacement field D, and the polarization P 
are well defined inside the slab. In practice, these fields 
may be calculated, for example, from unit cell avetages of 
the electrostatic potential and the charge density.B How- 
ever, we will see that for slabs with only 7-9 atomic lay- 
ers, the applicability of this approach has its limitations. 

In the case of an applied external electric field Ecxt 
perpendicular to the surfaces, the dielectric displacement 
field D inside the slab is oriented parallel to the z-axis 
and is equal to Eoxt- The boundary condition of a vanish- 
ing external electric field is therefore equivalent to a van- 
ishing dielectric displacement field D inside the slab. Fig- 
ure ||(a) shows a schematic picture of the planar-averaged 
potential v{z) for this situation. The potential is constant 
outside the slab, but due to the slab dipole moment m, 
the potential jumps by 47rem when going from one side of 
the slab to the other. At the same time, the polarization 
P leads to surface charges a—P ■ n which give rise to a 
huge depolarization field E=D— 47rP=— 47rP inside the 
slab (notice that E does not depend on the thickness of 







£ext=0 ^ 




A%e m 


depolarisation / 1 ^ 

field ---- _ / n 


^ 1^ £ext=0 


yv\j 


\lf- 


z 




FIG. 1. Schematic illustration of the planar-averaged po- 
tential v{z) for an isolated slab with a dipole moment m 
perpendicular to the surface, (a) Vanishing external electric 
field (equivalent to D—Q). (b) Vanishing internal electric field 

(£;=o). 



the slab). The contribution of the depolarization field to 
the total energy is Ipjg^ enough to completely destabilize 
the bulk FE state.ll3ll3 Therefore, relaxing a polarized 
slab under the boundary condition of a vanishing exter- 
nal electric field will inevitably result in a paraelectric 
cubic structure. 

A comparison with the situation in an infinitely ex- 
tended crystal shows which boundary conditions have 
to be used instead. In DFT calculations for bulk sys- 
tems, periodic boundary conditions for the electrostatic 
potential are usually applied. In this case, the internal 
electric field E vanishes, even in the presence of a sponta- 
neous polarization P^""', whereas the dielectric displace- 
ment field D will be nonzero. In analogy, a slab is in 
a FE state with spontaneous polarization Pg if a situa- 
tion exists where the internal electric field E is zero, but 
the dielectric displacement field D is nonzero. However, 
an external electric field Ecxt=D=47rPs will then appear 
outside the slab, as shown in Fig. ||(b). 

In our calculations we are only able to control the ex- 
ternal electric field Eoxt but not the internal field E. So 
to study whether a slab shows a FE instability, we have 
to apply external electric fields of different strength and 
search for the situation where E=Ecxt~47rP is zero. A 
rough estimate of how large Eoxt will be in this situa- 
tion can be made if we assume that the polarization of 
the slabs is equal to the bulk spontaneous polarization 




FIG. 2. Sketch of the Gibbs free energy G and the internal 
electric field i5 as a function of the dielectric displacement D 
(or the external electric field i?ext, respectively) for a para- 
electric (dotted line) or ferroelectric (solid line) slab. 

pbuik_ gy using the Berry phase approachJlZl we have cal- 
culated the bulk spontaneous polarization of BaTiOs and 
PbTiOa to be 4.6x10-^ e/bohr^ and 14.2x10-3 e/bohr^, 
respectively. This translates to external electric fields of 
0.058 a.u. ancLO.lS a.u. (atomic units are used through- 
out the paperta). 

To be able to distinguish between a paraelectric or fer- 
roelectric behavior of a slab, we employ a simple phe- 
nomenplpgical picture. We introduce the elastic Gibbs 
energyEj G(D) and we concentrate only on the depen- 
dence on the dielectric displacement field D (neglecting 
temperature and strain effects this is just the internal en- 
ergy U). For a paraelectric, G{T)) is roughly quadratic 
in D, whereas for a FE material GiTi) will show the well 
known double- well structure (see Fig. H(a)). Differenti- 
ating the Gibbs function gives immediately the internal 
electric field 



E 



dC 



(4) 



Thus, calculating E(D) directly reveals whether a slab is 
FE or not (see Fig. H(b)). For a paraelectric slab, E is 
proportional to D, with the proportionality factor given 
by the reciprocal dielectric constant, for a ferroelectric, 
the internal electric field is first negative as the polariza- 
tion builds up, and then vanishes when the spontaneous 
polarization Pg is reache d. W e will make use of this kind 
of analysis later in Sees. IIIB and [II C 3. 



2. Periodically repeated supercells 

In supercell calculations, periodic boundary conditions 
are usually imposed on the electrostatic potential. For 
slabs with a non- vanishing dipole moment perpendicular 
to the surface, this leads to electrostatic potentials that 
typically look like the sketch shown in Fig. ||(a). The 
electrostatic potential corresponds neither to the situa- 
tion of Fig. Wa) nor Fig. 0(b). Instead, the imposition of 
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FIG. 3. Schematic picture of the planar averaged potential 
v{z) for periodically repeated slabs, (a) with periodic bound- 
ary conditions, (b) potential of the dipole layer, (c) dipole cor- 
rected slabs with vanishing external electric field, (d) dipole 
corrected slabs with vanishing internal electric field. 



periodic boundary conditions on the supercell geometry 
leads to some other particular combination of internal 
and external electric fields, such that there is no discon- 
tinuity in the potential at the supercell boundary. The 
same occurs even for paraelectric slabs when terminated 
by nonequivalent surfaces with different work functions. 
The artificial electric fields become smaller when the 
thickness of the slab or the vacuum region is increased, 
but it is computationally very expensive to converge re- 
sults by using larger and larger supercells. Fortunately, 
the error associated with the artificial electric field can 
easily be eliminated by introducing an exteriial-dipole 
layer in the vacuum region of the supercell.ESEj The 
electrostatic potential of this dipole layer is shown in 
Fig. ^b). In order to reach a situation corresponding 
to Fig. JHa), the unwanted artificial external electric field 
can be compensated by adding a certain amount of the 
potential of Fig. p|(b) to that of Fig. p|(a), as is shown in 
Fig. |3|(c) (dipole correctionlla) . The dipole-corrected elec- 
trostatic potential is now discontinuous, but the discon- 
tinuity lies in the vacuum region of the supercell where 
the wave functions are essentially zero. Alternatively, the 
external dipole layer may be used to apply a true exter- 



nal electric field Ecxt to the surfaces£3 In particular, the 
situation of Fig. |l|(b) can be reached by subtracting a cer- 
tain amount of the dipole potential of Fig. ^(b) from that 
of Fig. ||(a), as shown in Fig. ^(d). 

The external dipole field can easily be implemented in 
any plane- wave based electronic, structure code. Follow- 
ing the notation of Bengtssonjlj we denote the external 
dipole potential of Fig. ^(b) as u'^'p(z), and the electro- 
static potential for the electrons calculated under peri- 
odic boundary conditions (corresponding to Fig. ^(a)) as 



'(r). The new potential is then 



(5) 



For a slab with dipole moment m and an external electric 
field -Eext, the dipole potential is given by 



„dip 



(^) 



( A'Km 

-e i/oxt 

V Co 



-y < ^ < y , (6) 



where cq is the height of the supercell. In a self-consistent 
calculation, the charge density, and thereby ttt., change 
with each step of the iteration. Therefore, m and the 
dipole potential v'^^^{z) have to be recalculated on each 
iteration until self-consistency is achieved (analogous to 
the updating of the Hartree and the exchange-correlation 
potentials) . 

The additional external potential w'^'p(z) also leads to 
changes in the total energies i?tot and the Hellmann- 
Feynman forces F/: 



E, 



tot — -Bjot 



27r77l 

Co 



-E,. 



A^ 



pper 



eZt 



f A'Km 
\ Co 



E, 



(7) 



(8) 



where S^T and Fj^"^ are the total energy and Hellmann- 
Feynman force calculated with the periodic potential 
tiP"'(r), and Z/ is the ionic charge of ion /. 



III. RESULTS AND DISCUSSION 

A. Zero external electric field 

As a first step, we calculated the fully relaxed structure 
of the various slabs in zero external electric field. For the 
asymmetrically terminated 8-layer slabs, the dipole cor- 
rection was used to enforce the vanishing of the e xterna l 
field. In this case, as has been pointed out in Sec. II C 1 , 
all slabs adopt the paraelectric cubic phase. The sym- 
metrically terminated slabs show no dipole moment (the 
central layer of the slabs is a mirror plane), and the small 
dipole moment of the asymmetrically terminated slabs is 
solely caused by the difference between the work func- 
tions of the two surfaces. 
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FIG. 4. Schematic illustration of the structure of the first 
three surface layers. 



Detailed results on the structure of the BaTiOs and 
PbTiOa surfaces in the cubic phase in the absence of 
an external electric field have already been published in 
Refs. S and 0. Therefore, we give here only a brief sum- 
mary in order to establish notation and to provide a base- 
line for comparison. 

Figure ^ shows a schematic illustration of the structure 
of the first three surface layers, where the relaxations 
of the atoms have been highly exaggerated. For both 
BaTiOa and PbTiOa, and for both the BaO/PbO and 
Ti02 surface terminations, the same characteristic fea- 
tures appear regarding the buckling of the surface lay- 
ers and the changes of the interlayer distances. Only 
the amplitudes of these relaxations differ from surface to 
surface. For a quantitative analysis of the surface relax- 
ations, we let 5z{0i) and 5z(yi-i) be the displacements 
of the oxygen and metal atoms, respectively, relative to 
the ideal unrelaxed structure in layer i. The change of 
the interlayer distance Ady is then given by the differ- 
ence of the averaged displacements {5 ziyi-)~\-5 z{0)) / 2 of 
the atoms in layers i and j. For all surfaces, A(ii2 is 
negative (corresponding to a reduction of the interlayer 
distance between the first and the second surface layer 
compared to the bulk value), whereas Ad23 is positive. 
To describe the buckling of the surface layers, we define 
a rumpling parameter rji as the amplitude of the relative 
displacements between the metal and the oxygen ions: 
{5ziyi-i)—5z{Oi))- "Hi is negative if the metal ions are be- 
low the oxygen atoms, which is true for the first surface 
layer of all surfaces. For the next surface layers, rji os- 
cillates in sign from layer to layer, and the amplitude 
decreases very rapidly. 



B. External electric field without field-induced 
atomic relaxations 

In the next step, we exposed the slabs to external elec- 
tric fields Ecxt of increasing strength, but we kept all 
atoms frozen in the positions which have been calculated 
in zero external field. Only the distribution of the elec- 




BaTiO, 



FIG. 5. Planar-averaged screening charge distribution of 
a BaO-terminated slab of BaTiOa , calculated from a 7-layer 
slab for an electric field of _Ecxt=+0.02 a.u. The arrows indi- 
cate the positions of the atomic planes. "Top" of slab is at 
right. 

trons was recalculated self-consistently. This mimics the 
exposure of the slabs to a AC electric field with a fre- 
quency high enough that the ions are no longer able to 
follow. In this situation, we expect the slabs to show a 
purely paraelectric response. That is, a linear relation 
between the internal electric field E and the dielectric 
displacement field D should be found. D=eooE, where 
Eoo is the optical dielectric constant. 

To verify this, we calculated the internal electric field 
E as a function of the applied external electric field Eoxt 
(or, equivalently, the dielectric displacement field D). A 
natural way to compute E would be to determine the gra- 
dient of the macroscopically averaged electrostatic poten- 
tial v{r). Unfortunaly, it turned out that 7-9 layers are 
not enough to give accurate macroscopic averages (the 
uncertainties in the gradients were much too large). In- 
stead, we calculated the internal electric field by using 
the relation E=D— 47rP. The polarization P can either 
be deduced from the surface charge a or from the dipole 
moment m and the thickness d of the slabs: 



P = (J , 



P = m/d 



(9) 



The dipole moment m is directly given by the charge dis- 
tribution of the slabs via Eq. (||) . For the thickness d we 
have taken the distance between the centers of gravity zq 
of the screening charge distribution of the top and bottom 
surface of our slabs. (For a metal, zq is the position of the 
surface image plane from which the classical image po- 
tential is measured.) The screening charge distribution is 
given by the difference of the charge densities calculated 
with and without an external electric field. A typical 
planar-averaged screening charge distribution p^'^''[z) is 
shown in Fig. |^. As expected, a positively-oriented elec- 
tric field pushes the electrons into the slab at the top sur- 
face and pulls them out at the bottom (where it can be 
interpreted as a negatively-oriented field). From Fig. H, 
we see that the screening charge piles up directly at the 
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FIG. 6. Internal electric field E in the 7-layer slab with 
frozen atomic positions as a function of the applied external 
electric field -Ecxt. Solid lines: AO-terminated slabs, plotted 
lines: Ti02-terminated slabs. Atomic units are usedJU 
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TABLE I. Calculated optical dielectric constants Eoo for 
the various slabs. The theoretical bulk values are taken from 
Ref. Eil. 



surface and is mainly confined above the surface atoms. 
These surface charges are responsible for the screening 
of the external electric field. Inside the slab, fairly large 
oscillations are found which are associated with the re- 
maining internal electric field. 

The surface charge a and the center of gravity of the 
screening charge distribution zq for the top surface of the 
slab are given by 



co/2 
cr = / p""(z) dz , 

Zcut 



co/2 

zn = ^ / p^"(z)zdz . (10) 



For Zcut we have taken the position where the extrapola- 
tion of the first charge peak goes to zero (see Fig. ||), but 
the results are not very sensitive to the choice of Zcut- 
Analogous relations apply for the bottom surface. To 
check the consistency of this approach, we have calcu- 
lated the polarization P of the slabs in both of the ways 
indicated in Eq. (0), i.e., using either the surface charges 
a or the dipole moment m with the thickness d. In all 
cases, the computed values for the polarization differ by 
less than 1%. 

With the polarization P we can finally determine the 
internal electric field E. Figure ^ shows the internal elec- 
tric field for the 7-layer slabs as a function of the applied 
external field. Compared with our phenomenological pic- 
ture in Fig. 0, a clear paraelectric behavior can be seen 
for all slabs. Table | gives the optical dielectric constants 
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FIG. 7. Field dependence of the rumpling parameter r\i calculated from the 7-layer slab. Dotted horizontal lines: am- 
plitude of displacement of metal and oxygen ions in the bulk ferroelectric state. Dotted vertical lines: External electric field 
corresponding to the bulk spontaneous polarization. 



deduced from the E{U) plots for the different slabs, to- 
gether with the LDA bulk values of Ref. |l]. Already 
for 7 atomic layers, the dielectric constants agree very 
well with the bulk values. The finite-size effect is less 
than 10%. Its sign depends on the relative numbers of 
AO and Ti02 layers: Coo is increased when the number 
of AO-layers exceeds the number of Ti02-layers, and is 
reduced when the slabs contain more Ti02 layers. 



C. External electric field with atomic relaxations 

Finally, all slabs were fully relaxed in the presence of 
external static electric fields. Fields of strength varying 
between —0.08 a.u. and -t-0.08 a.u. have been applied. 
Starting from zero field, the amplitude of the field was 
changed in steps of 0.01 a.u., and the atomic positions at 
the previous field were used as initial configurations for 
the new atomic relaxations. The upper limit of 0.08 a.u. 
for the amplitude of the external electric field is deter- 
mined in our calculations by the thickness of the vac- 
uum region for the following reason. As can be seen in 
Fig. P(d), there is a kind of "quantum well" located in 
the vacuum region just to the right of the dipole layer. 
If the electrostatic potential of this quantum well drops 
below the Fermi level, it can become populated by the 
transfer of electrons from the slab region. The threshold 
for the occurrence of this unwanted behavior depends on 
both the width of the vacuum region and the strength of 
the electric field. 



1. Field-induced atomic relaxations 

In a positively oriented external electric field, the neg- 
atively charged oxygen ions will be pushed into the bulk 
and the positive ions will be pulled towards the surface. 
This will immediately change the rumpling parameter rji 
of the surface layers, as can be seen in Fig. 0. For small 
applied fields, rji increases linearly with the strength of 
£^oxt, as is expected if a harmonic coupling between the 
atoms is assumed. However, for larger fields, nonlinear 
effects become important, especially for the first and sec- 
ond AO layers. 

To facilitate comparison of these surface rumplings 
with the corresponding bulk rumplings of the "up" and 
"down" FE states, the latter are drawn as horizontal dot- 
ted lines in Fig. 0. That is, these lines show the values 
that the rumpling parameters would have if the corre- 
sponding slab had been a piece of truncated FE bulk 
material without any further relaxation of the atoms. 
Additionally, for BaTiOs we have used dotted vertical 
lines to indicate the magnitude of the external electric 
field Eext=47rPg""^ that would be consistent with a bulk 
spontaneous polarization P^""^ in the interior of the slab. 

It can be seen from Fig. |^ that for this external electric 
field, the rumpling parameters r]i of the surface layers 
have reached the same order of magnitude as the bulk 
FE distortions. For most layers rji is slightly larger, but 
for some layers it is a little bit smaller. Overall, the 
rumpling parameter rji tends to be larger for the nega- 
tive electric field than for the positive field, suggesting 
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FIG. 8. Field-dependence of the change in interlayer dis- 
tances Adij calculated for the 7-layer slabs. 



that the surface polarization is larger for the negatively- 
oriented field. However, it is difficult to deduce from 
Fig. R whether the polarization is enhanced or reduced 
at the surface compared with the bulk. We will return 
to this question in the next subsection. 

The bulk spontaneous polarization P^""^ is roughly 
three times larger for PbTiOa than for BaTiOa (see Sec- 
tion If C f ) . The external electric field for the case that 
a slab adopts the PbTiOs bulk spontaneous polarization 
is therefore on the order of 0.18 a.u. This is well beyond 



the upper limit of the electric field that can be applied in 
our calculations. Consequently, the polarization of the 
PbTiOs slabs is well below P^""^ for all applied exter- 
nal electric fields, and therefore it is no surprise that the 
surface rumpling iji in Fig. is much smaller than the 
corresponding bulk FE rumpling over the entire range of 
accessible external electric fields. 

Associated with the changes in the layer rumplings, we 
also observe changes in the interlayer distances as shown 
in Fig. p. (Note that the atomic layers are charge-neutral, 
so that there is no net force on the layers in an external 
electric field.) We find some noteworthy qualitative dif- 
ferences in the behaviors of the two materials (e.g., in 
the slopes of the Ad23 vs. E'cxt curves). Ultimately these 
differences can presumably be traced to the rather differ- 
ent chemistry of Ba-0 and Pb-0 bonds, but a detailed 
explanation of the observed trends is not obvious. 
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FIG. 9. Dipole moment per atomic layer calculated for 
the 7-layer slabs for electric fields from to 0.08 a.u. in steps 
of 0.01 a.u. The orientation of the electric field corresponds to 
a positive field for the top surface and a negative field for the 
bottom surface. The dotted line represents the dipole moment 
for an atomic layer in the bulk ferroelectric state, calculated 
from the bulk spontaneous polarization via c/2-P^^^^. 



2. Field-induced layer dipoles 

To address the question how the polarization is 
changed at the surface, we have subdivided the total 
dipole moment m of the slabs from Eq. (m) into dipole 
moments 771^ per atomic layer. This was done by iden- 
tifying the "nodes" Zi at which J' ,„ p{z) dz vanishes. 
Regions between successive nodes are thus charge neutral 
and contain exactly one AO or Ti02 atomic plane. The 
dipole moment rrii was then calculated via Eq. dl^) while 
restricting the range of integration to the area between 
the two nodes which embrace the atomic layer of interest. 

Fig. shows the dipole moments per atomic layer for 
the 7-layer slab in different external electric fields. For 
zero external electric field, the dipole moments oscillate 
in sign from layer to layer. The dipole moment is largest 
for the first atomic layer at the surface and goes to zero 
very rapidly (see also Fig. Q). Finally, the net dipole 
moment of the surface, given by the sum of the mi of 
all surface layers, points inwards towards the bulk for all 
surfaces. 

An external electric field causes the dipole moments of 
all layers in the slab to increase more or less equally. The 
dipole moments continue to oscillate from layer to layer 



dipole moment 
AO-term. Ti02-term. 



polarization 
AO-term. Ti02-term. 



BaTiO, 



PbTiO, 



SrTiOa -8.7 -14.2 1.19 (-) 1.95 (-) 

BaTiOa -7.4 -14.7 0.98 (21%) 1.95 (42%) 

PbTiOs -11.3 -11.8 1.48 (10%) 1.55 (11%) 

TABLE II. Dipole moment (in 10~^e/bolir) and polariza- 
tion (in 10~"^e/bohr^) pof the surface unit cell as calculated 
from the 7-layer slabs.tj The surface polarization as a frac- 
tion of the theoretical bulk spontaneous polarization is given 
in parentheses. Values for SrTiOa are given for comparison. 



when i^oxt is increased, but they approach a constant 
bulk value already after two layers in BaTiOs and after 
three layers in PbTiOs. (This becomes more apparent 
from the 8- and 9-layer slab calculations, not shown here.) 
The dipole moments in the surface layers arc modified 
relative to the bulk layers by the inwards oriented surface 
dipoles. This leads to a suppression of the polarization at 
the surface in a positively-oriented external electric field, 
and an enhancement for a negative field. However, the 
modification of the polarization is more or less confined to 
the first unit cell (i.e., first pair of layers) at the surface. 

In most cases a rough estimate of how much the polar- 
ization is suppressed or enhanced at the surface is given 
by the surface dipole moments in zero external electric 
field. Table |l| lists the dipole moment of the uppermost 
unit cell (i.e., the sum of the dipole moments rrii of the 
first two atomic layers) in zero external electric field for 
slabs of 7 layers. (The results from 8- and 9-layer slabs 
differ by less than 3%.) Additionally, these dipole mo- 
ments have been converted to a polarization by dividing 
them by the lattice constant c. 

The values of the dipole moments of the uppermost 
unit cells are of the same order of magnitude for all sur- 
faces. In the case of BaTiOa, they correspond to about 
20-40% of the bulk spontaneous polarization Pg"'"*^. 
Therefore, for BaTiOs a relatively large change in the 
polarization at the surface has to be expected for a spon- 
taneously polarized slab. On the other hand, for PbTiOs 
the surface dipoles make up only about 10% of the bulk 
spontaneous polarization. The modification of the polar- 
ization at the surface in a spontaneously polarized state 
will therefore be only modest. The polarization will be 
suppressed for an upwards polarization and enhanced for 
a downwards polarized state. 



3. Ferroelectric instability of the slabs 

The results presented so far on the structure of the sur- 
faces do not provide an answer to the question whether 
the slabs exhibit a true FE instability or simply show 
a paraelectric polarization. Instead, we have to deter- 
mine the electric field inside the slabs as a function of 
the applied external field and compare the result with 
our simple phenomenological picture from Fig. 0(b). 



o 

til 




-1.0 



-0.08 -0.04 0.00 0.04 0.08 -0.08 -0.04 0.00 0.04 0.08 

^ax, [a-U-1 £ex, [a-u] 



FIG. 10. Internal electric field E in the fully relaxed 
7-layer slab as a function of the applied external electric 
field -Eoxt. Solid lines, AO-terminated slabs; dotted lines, 
Ti02-terminated slabs. 



For the calculation of the internal electric field wc fol- 



lowed the same procedure as described in Section IIIB. 
However, instead of using the true screening charge dis- 
tribution (which includes contributions from the moving 
ions) to determine zq and d, we took the difference be- 
tween the charge density of the fully relaxed slab in the 
external electric field and the charge density of a calcu- 
lation where the atomic positions have been kept fixed 
and the electric field turned off. Unfortunately it turned 
out that it is not possible to determine the internal elec- 
tric field to very high precision. The use of the relation 
E=D— 47rP involves the difference of two relatively large 
numbers, so that the result for the internal electric field 
depends very sensitive on details of the calculation. In 
particular, the determination of the thickness d of the 
slabs is very critical. Also when approaching the break- 
down field above which electrons accumulate in the quan- 
tum well next to the dipole layer, the wavefunctions start 
to penetrate the barrier between the slab and the quan- 
tum well. This leads to a small overestimation of the 
dipole moment and the polarization of the slab at high 
applied fields. But because of the subtle difference be- 
tween the external electric field and the polarization, the 
error in the internal electric field becomes noticable, and 
deviations from the behavior described in Fig. ^(b) may 
appear at large external fields. Therefore, the curves of 
the internal electric field as a function of the applied ex- 
ternal electric field shown in Fig. Uu represent only basic 
tendencies but are not to be taken as quantitatively ac- 
curate results. 

The results in Fig. 10 are clearest for the BaO- 
terminated slab of BaTiOs, where the full curve ex- 
pected for ferroelectric behavior, as in Fig. 0(b), can 
be observed. We find that the internal electric field 
vanishes for an external electric field of approximately 
0.05 a. u. (the result for the 9-layer slab is the same). This 
translates to a spontaneous polarization of the slab of 
Ps=4x 10"'^ e/bohr^, which corresponds roughly with the 
bulk spontaneous polarization P^""^ of BaTiOs. How- 



ever, due to the limitations in the accuracy of our cal- 
culations, we are not able to say for certain whether the 
spontaneous polarization of the slab is enhanced or re- 
duced compared to the bulk value. 

For the other three cases in Fig. nO it is more difficult 
to deduce, whether the slabs show a FE instability or not. 
For PbTiOs we have the problem that the upper limit for 
the external electric field is reached long before the slabs 
come close to a polarization equal to the bulk sponta- 
neous polarization (as has been pointed out above). In 
Fig. O, the external electric field where the internal field 
is zero can therefore not be reached for the PbTiOa slabs. 

On the other hand, according to Fig. p|(b), a nega- 
tive slope of the internal electric field at zero applied 
field is already a very good indicator for a FE instability. 
As can be seen from Fig. RQ this condition is fulfilled 
for both the PbO-terminated and the Ti02-terminated 
slab of PbTiOa, but not for the Ti02-terminated slab 
of BaTiOa. To estimate how sensitive this result de- 
pends on details of our calculations we compare the nu- 
merical values of the initial slopes of the E{D) curves 
with an error bar which we derive from the uncertain- 
ties in the determination of the slab thickness d. The 
fit curves in Fig. Il^ yield initial slops of —9.1x10"^ and 
-f 2. 6x10^3 foj. tjjg Ba,0- and Ti02-terminated BaTiOs 
slabs and -2.5x10^3 ^^d -7.6 xlO'^ for the PbO- and 
Ti02-tcrminatcd PbTiOa slabs, respectively. By vary- 
ing Zcut in Eq- (|lO|) we can estimate an error bar of about 
±4x 10^^. For the BaO-terminated BaTiOs slab and the 
Ti02-terminated PbTiOs slab the initial slopes are well 
outside the error bar, indicating a ferroelectric instabil- 
ity of the slabs, whereas for the Ti02-terminated BaTiOa 
slab and the PbO-terminated PbTiOs slab our analysis 
of the internal electric field is not accurate enough to 
make a clear statement. It would be possible to clarify 
the situation in the last two cases by going to thicker 
slabs with more atomic layers, since then the thickness d 
would be relatively less uncertain, but this was deemed 
too computationally demanding to be undertaken in the 
current project. At least we can say that both slabs are 
very close to a ferroelectric instability. 



IV. SUMMARY 

By using first-principles density-functional calcula- 
tions, we have studied (OOl)-oriented slabs of BaTiOa 
and PbTiOa in external electric fields. Artificial elec- 
tric fields introduced by the periodic boundary conditions 
of the supercell approach have been eliminated with the 
help of an external dipole layer in the vacuum region. In 
this way it is possible to handle periodically repeated po- 
larized slabs as if they were genuinely isolated, without 
the need for large vacuum separations. This makes the 
dipole correction a very useful tool for all problems where 
slabs with nonvanishing dipole moments or inequivalent 
surface terminations have to be considered. 



For the BaO-terminated BaTiOa slabs and the Ti02- 
terminated PbTiOs slabs as thin as 7 atomic layers, we 
have found strong evidence for the presence of an in- 
stability to a spontaneously polarized ground state with 
the polarization perpendicular to the surface. This is in 
agreement with recent experimentsa and a microscopic 
effective-Hamiltonian studyJZl However, because the lat- 
ter approach does not explicitly take into account the 
structural relaxations at the surface, it cannot provide 
the detailed level of description that our theory gives. In 
particular, we find an enhancement of the polarization 
at the surface for which the polarization points inwards, 
whereas the polarization is reduced when the sponta- 
neous polarization points outwards. We have also de- 
tailed the variations in these surface relaxations as the 
polarization of the slab is modified by an external elec- 
tric field. 

In order to obtain an even more realistic level of de- 
scription, the next obvious step will be to calculate the in- 
fiuence of a real metal/perovskite interface on the atomic 
relaxations and the polarization at the interface. How- 
ever, such a study requires the choice of a particular 
metallic overlayer with particular epitaxy conditions, and 
is thus intrinsically less universal. Nevertheless, it may 
be an interesting avenue for future investigation. 
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